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ABSTRACT 

The observed characteristics of the nebulosity surrounding the SMC High Mass X- 
ray Binary XTE J0111. 2-7317 are examined in the context of three possible nebular 
types: SNR, bowshock and HII region. Observational evidence is presented which 
appears to support the interpretation that the nebulosity surrounding XTE J0111.2- 
7317 is an HII region. The source therefore appears to be a normal SMC Be X-ray 
binary (BeXRB) embedded in a locally enhanced ISM which it has photoionised to 
create an HII region. This is supported by observations of the X-ray outburst seen 
with BATSE and RXTE in 1998-1999. It exhibited characteristics typical of a giant 
or type II outburst in a BeXRB including large spin-up rates, Lx > 10 38 erg/sq.cm-s, 
and a correlation between spin-up rate and pulsed flux. However, the temporal profile 
of the outburst was unusual, consisting of two similar intensity peaks, with the first 
peak of shorter duration than the second. 



1 INTRODUCTION AND BACKGROUND 

The X-ray transient XTE J0111. 2-7317 was first detected 
by the Proportional Counter Array (PCA) , on the Rossi X- 
ray Timing Explorer (RXTE) X-ray observatory, on October 
the 29th 1998 (Chakrabarty et al 1998a; Chakrabarty et al 
1998b) in the 2-10keV band. 

The observed temporal, flux and spectral characteristics 
are typical of BeXRBs and the detection of pulsations at 31 
seconds further strengthens this hypothesis. In addition the 
source location is coincident with the SMC, the assumption 
of membership permitting an absolute luminosity to be cal- 
culated. Whilst towards the upper end of the typical flux 
distribution, peaking at approximately 2 x 10 38 ergs' 1 , it is 
consistent with a giant (Type II) outburst from a BeXRB. 

The optical counterpart was identified by Israel et al 
(1999) and subsequently optical and IR measurements by 
Coe, Haigh and Reig 2000 (henceforth referred to as CHR) 
confirmed the optical counterpart to be a B0-B2 luminos- 
ity class III-V star showing a strong IR excess. Surpris- 
ingly, however, their Ha images of the field revealed an 
extended region of Ha emission surrounding the source. 
They suggested it could be a surrounding SNR or wind bow 
shock. Otherwise XTE J0111. 2-7317 is in many ways typ- 
ical of the growing population of known SMC BeXRBs. If 
a SNR the size of 3.1pc, using the diameter-to-age relation- 
ship D = 0.9t 2/s (Clark & Caswell 1976), implies an age of 
120yrs. A typical expansion velocity of lOOOOkm/s (Cheva- 
lier 1977) for Type II SNe implies an age of ~ 300years. 
Such a young BeXRB would be of particular interest. 



Throughout this work XTE J0111. 2-7317 refers both to 
the X-ray source and to the optical counterpart. 



2 OBSERVATIONS 
2.1 Optical images 

Initial photometry was reported by CHR who carried out 
Johnson, Stromgren and IR measurements. For this work 
further Wide Field Imager (WFI) images of the field taken 
on 26 July 1999 were extracted from the ESO archive. The 
WFI is a focal reducer-type camera which is permanently 
mounted at the Cassegrain focus of the 2.2-m MPG/ESO 
telescope at La Silla, Chile. The detector is a mosaic of 
4 x 2 4096 x 8192 CCDs with 0.24" pixels. These obser- 
vations targeted a different object, but because of the large 
34' x33' field XTE J0111.2-7317 was also imaged. Useful im- 
ages in [SII], Ha, [OIII] and continuum regions close to 
Ha were de-archived. Preprocessing was performed using 
STARLINK software. 



2.2 IR measurements 

The source was observed as part of the DENIS IR survey 
of the Magellanic Clouds in Sept 1996. This IR catalogue 
shows I=15.27±0.02 and J=15.19±0.10 measurements, but 
not Ks, because of the source's faintness. However, a Service 
Observation of the field was carried out using IRIS2 on the 
AAT during commisioning time. The resulting K band image 
was calibrated using three DENIS sources in the field and 
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Figure 1. Ha region of SAAO XTE J0111.2-7317 spectrum with 
2 component Gaussian fits. 



an accurate determination of K=14.79±0.03 was achieved. 
The resulting J-K value of 0.4 is typical of other BeXRB 
systems in the SMC and indicative of an excess IR flux over 
the stellar continuum for the B0-B2 III-V star classified by 
CHR. 



2.3 Optical spectroscopy 

A blue spectrum was obtained at the ESO 1.52m telescope 
at La Silla, Chile, on November the 3rd 1999. The Boiler 
and Chivens spectrograph was used with the no: 33 holo- 
graphic grating, covering the range from the classification 
region (4000 - 5000A) up to ~ 5900A with a resolution of 
lA per pixel. The slit was oriented in a NW-SE line to en- 
compass the brightest parts of the nebulosity. 

CHR reported results from spectroscopy in which the 
only features seen were strong Ha and Hj3 emission lines, 
as expected for a BeXRB. The measured equivalent widths 
were -27±0.3A for Ha and -3.8±0.2A for Hf3. However, 
a considerable fraction of this emission is contamination 
from the nebular emission which crosses the stellar spec- 
trum. Because the nebulosity is brightest close to the star, 
the sky /background subtraction carried out as a standard 
procedure in the spectral reduction does not remove such 
contamination. The nebular [SII] lines at 6716A and 673lA 
contaminate the stellar spectrum in the same way. All pub- 
lished equivalent widths to date (CHR and Covino et al. 
2001) have been skewed by this effect. 

Figure 1 shows the Ha region of the SAAO XTE 
JO 11 1.2- 73 17 spectrum reported in CHR. Inspection of the 
spectral images shows that while the nebular contamina- 
tion is fully contained in a band 5 pixels wide (~ 2.1 A 
or 97km/s), the true circumstellar emission is considerably 
broader at ~ 400km/s, typical for Be stars. Thus whilst the 
measured EWis certainly an overestimate, there is no doubt 
that intrinsic circumstellar emission exists. A two compo- 
nent Gaussian fit to the line profile was performed within 
the DIPSO/ELF package using the measured FWHM of the 
nebular contamination (obtained at a spatial position adja- 
cent to the star) to robustly remove it (Figure 1). This fit 
assigns 61 ± 4% of the line flux to the broader circumstel- 
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Figure 2. The 4800-5050A region of the ESO spectrum of XTE 
J0111. 2-7317 showing H/3, [OIII] and Hel emission lines. The 
same range of a spectrum of the classic BeXRB X Per from the 
same observing run is shown for comparison. The XTE JO 11 1.2- 
7317 spectrum has been shifted by — 165km/s to the heliocentric 
rest frame. 



lar disc (conventional Be) component, producing an intrinsic 
EWHa of 16.5±l.lA. 

The Doppler shift of the broad circumstellar component 
of Ha is recessional at 170.8±9.3km/s (visible in Figure 1) 
which agrees well with the systemic value of 166±3km/s 
(Feast, 1961) for the SMC. Not only does this confirm mem- 
bership of the SMC (never really in doubt because of the 
X-ray flux and photometry) but it also places a constraint 
on the system's radial velocity relative to the local ISM of 
2.5 ± 12.5km/s. Accurate Doppler shifts for all of the lines 
detected from the surrounding nebulosity are in Table 2. 

Regarding a classification based upon spectroscopic fea- 
tures, this has been done by Covino et a. (2001) based upon 
superior spectra, and their conclusion of B0.5-lVe is consis- 
tent with this work. 

Spectra (Figure 2) show lines at 4959A and 5007A, at- 
tributed to [OIII] emission. Such forbidden line emission is 
typical of many nebulae but uncharacteristic of the rela- 
tively dense circumstellar environment of Be stars. A search 
through the spectra of 11 other LMC/SMC BeXRB spectra 
in the Southampton database, and several galactic systems 
(including X Per and A0535+262) found no other systems 
with features at this wavelength (as expected if it arises from 
[OIII] emission). It is shown below that this feature results 
from the superposition of nebular emission. 



2.4 X-ray data 

2.4-.1 Observations 

XTE J0111. 2-7317 was active in X-rays from 1998 Octo- 
ber - 1999 February. Figure 3 shows the 2-10 keV flux his- 
tory measured with the All-Sky Monitor (ASM) (Levine 
et al, 1996) on the RXTE. XTE J0111.2-7317's location 
was in the field-of-view of the RXTE PCA (Jahoda et 
al, 1996) and High Energy Timing Experiment (HEXTE) 
(Rothschild et al.1998) for 42 observations from 1998 Oc- 
tober 20 - 1999 May 11. Further, the Burst and Transient 
Source Experiment (BATSE) (Fishman et al, 1989) on the 
Compton Gamma Ray Observatory had also observed XTE 
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Figure 3. 2-10 keV flux history for XTE J0111.2-7317 measured 
with the RXTE ASM. All points are 4-day averages. 



J0111.2-7317's location from 1991-2000. No additional out- 
bursts have been detected with either BATSE or RXTE. 



Figure 4. Top: Barycentered pulse frequencies measured with 
BATSE (filled squares) at 4-day intervals and measured with the 
RXTE PCA for each observation (open circles). Center: 20-50 kcV 
pulsed flux measured with BATSE at 4-day intervals. (1 mCrab = 
1 X 10 -11 ergs cm -2 s _1 ) Bottom: 2-60 kcV pulsed flux measured 
with the RXTE PCA for each observation. (1 mCrab = 2.786 cts 
s- 1 PCIJ- 1 ) 



2.4.2 Timing 

To determine pulse frequencies for XTE J011 1.2-7317, we 
performed a grid search over a range of candidate frequencies 
using 1-second resolution 20-50 keV data from the BATSE 
Large Area Detectors (LADs). This technique is described 
in detail elsewhere (Finger et al, 1999; Wilson- Hodge 1999; 
Wilson et al, 2002; Wilson et al, 2003). First we combined 
the count rates over the 4 LADs viewing XTE J0111. 2-7317, 
using weights optimized for an exponential energy spectrum 
f(E) = Aex.p(-E/kT) with temperature kT = 12 keV, 
and then grouped them into 300-s segments. In each 300- 
s segment, we fitted a model consisting of a sixth-order 
Fourier expansion in pulse phase model (representing the 20- 
50 keV pulse profile) and a spline function with quadratics 
in time (representing the background) . Using 300-s segments 
allowed use of a quadratic background model, allowed us to 
effectively fit bright Earth occultation steps in the data, im- 
proved computational efficiency, and reduced smearing of 
the pulse profiles if the initial phase model was incorrect. 
The final step was to perform a grid search in frequency 
using the set of typically several hundred estimated 20-50 
keV pulse profiles in each 4-day interval of data. For each 
grid point, the segment pulse profiles were shifted accord- 
ing the the corresponding frequency offset. These profiles 
were then combined into a mean profile. The best frequency 
was selected using a modified Z% statistic, called Y n af- 
ter (Finger et al, 1999), given by Y n = Y2=i Ml<& h 
where olk is the mean Fourier coefficient for harmonic h and 
a% h is the sample variance from the segment profiles. Us- 
ing sample variances from the segment profiles rather than 
assuming Poisson statistics properly accounted for aperi- 
odic noise in XTE J0111. 2-7317 and in any other sources 
in the large BATSE field of view. A similar technique was 
used to generate pulse frequency measurements for RXTE 
PCA observations from 2-60 keV Standard 1 (125 ms, no 
energy resolution) data and from RXTE HEXTE data. The 
root-mean-squared (rms) pulsed fluxes were estimated at the 
best fit frequency as Frms = (0.5 27JhLi l^h | 2 ) 1//2 . Figure 4 
shows the barycentered pulse frequency history measured 



with BATSE (squares) and with the RXTE PCA (open cir- 
cles), the 20-50 keV pulsed flux measured with BATSE, and 
the 2-60 keV pulsed flux measured with the RXTE PCA. 
Pulsations were detected with BATSE from 1998 October 
30 - 1999 January 9. The RXTE PCA continued to detect 
pulsations until observations ceased on 1999 February 19. 
Pulsations were not detected in two subsequent PCA obser- 
vations on 1999 March 26 and 1999 May 11. 

Pulse frequency derivatives were computed by differ- 
encing adjacent BATSE pulse frequency measurements. Fig- 
ure 5 shows the correlation between frequency derivative and 
20-50 keV pulsed flux. This correlation is well fitted with a 
power law with index 0.9 ± 0.1, which is consistent with the 
index of 6/7 expected from simple accretion theory, if accre- 
tion from a disk is assumed. If a reliable bolometric correc- 
tion can be derived for a source, this correlation can be used 
to constrain the distance or the magnetic field strength if the 
other is known. Unfortunately in this case, large changes in 
the pulse fraction over the course of the outburst prevented 
us from deriving a reliable bolometric flux (see Figure 6). 



2.4-3 Energy Spectra 

Analysis of the energy spectrum of XTE J0111. 2-7317 was 
complicated by the fact that SMC X-l was also present in 
the RXTE PCA and HEXTE field-of-view in many of our 
observations. Using power spectra from each observation, we 
determined in which observations 31-s pulsations from XTE 
J0111. 2-7317 were present and 0.7-s pulsations from SMC 
X-l were not present. We then generated energy spectra and 
response matrices for the PCA and HEXTE data for these 20 
observations and fitted them in XSPEC with an absorbed 
power-law with a high energy cut-off and a Gaussian iron 
line. Average parameter values for the 20 observations are 
given in Table 1. The flux in the iron line was correlated 
with the total flux, with a correlation coefficient of 0.97 and 
a chance probability of 2.5 x 10~ 5 , suggesting that the iron 
line is intrinsic to XTE J0111. 2-7317. No other parameters 
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Figure 5. Frequency derivatives versus 20-50 keV pulsed flux 
measured with BATSE. The solid line denotes a power law fitted 
to these data with an index of 0.9 ± O.f . 

were obviously correlated with flux nor did they show clear 
evolution with time. 

Using only those observations where SMC X-l was not 
present in the power spectrum, we estimated pulse fractions 
from RXTE PCA data. Figure 6 shows the pulse fraction 
versus the 2-50 keV flux computed from our XSPEC fits. 
The pulse fraction is correlated with the total flux, indi- 
cating that the increase in flux is primarily in the pulsed 
component. 



3 DISCUSSION 

3.1 The X-ray observations 

Comparing BATSE and PCA pulsed fluxes, we see evidence 
that the first peak of the outburst is harder than the second 
peak, shown in Figure 7. Unfortunately, because SMC X-l 
is present in all of the observations during the first peak of 
the outburst, we cannot study the effects of this change in 
hardness on the shape of the energy spectrum in detail. The 
physical explanation of this is unclear. We speculate that 
this may be an obscuration or absorption effect - the softer 
X-rays are partially absorbed or blocked by the large amount 
of material present at the beginning of the outburst. As the 
outburst progresses, this material begins to dissipate and 
the softer X-rays become more easily seen. A comparison of 
HEXTE and PCA pulsed fluxes suggests a similar change in 
hardness, although the HEXTE fluxes are much less certain 
than the BATSE fluxes due to the smaller effective area and 
integration time. This suggests that the effect is intrinsic to 
the source and is most likely not instrumental. 

Observations of steady spin-up during the outburst al- 
low us to compute a lower limit to the luminosity of XTE 
J0111. 2-7317. The angular momentum of a rotating neu- 
tron star is given by I = 2-kIv where v is the spin fre- 
quency and / is the moment of inertia of the neutron star. 
The torque on the neutron star is obtained by differenti- 
ating £ or N = £ = 2-kIv. In an accreting pulsar under- 
going steady spin-up, the characteristic torque is given by 
iVmax = M(GMxrm) 1 ' 2 where M is the mass accretion rate, 
G is the gravitational constant, Mx is the mass of the neu- 
tron star, and r m is the magnetospheric radius. The max- 



Table 1. X-ray Spectral Parameters for XTE J0111. 2-7317 



Parameter 


Value 


Standard Deviation 


N H (cm" 2 ) 


1.5 x 10 22 


0.5 x 10 22 


index 


0.93 


0.06 


£ cut (keV) 


13.4 


0.8 


Bfoid (keV) 


20.3 


2.2 


E linc (keV) 


6.5 


0.2 


o"iinc (keV) 


0.8 


0.5 



imum possible torque occurs when the magnetospheric ra- 
dius equals the corotation radius. Setting N < .iV max yields 
an expression for M and hence Lx that depends only on the 
spin frequency and its derivative for assumed values of the 
neutron star parameters, i.e. 

Lx > (27r) 4/3 J(GM x ) 1/3 i? x 1 ^ 1/3 . (1) 

Assuming I = 10 45 g cm 2 , Mx = 1.4M , and Rx = 10 
km, Lx > 10 38 ergs s _1 for typical observed values of v — 
32.4 mllz and v = 5 x 10~ n Hz s _1 . The average 20-50 
keV pulsed flux was 2.5 x 10~ 10 ergs cm~ 2 s _1 , implying 
that XTE J0111. 2-7317 is at a distance of d > 59.8 kpc, 
confirming that it is in the SMC. 

The outburst from XTE J011 1.2-7317 should be classi- 
fied as a giant or Type II outburst given the large spin-up 
rates (3 - 7) x 10~ n Hz s _1 measured with BATSE and the 
large luminosities (2 — 9) x 10 38 ergs s _1 , measured in the 
PCA observations where SMC X-l was not present. How- 
ever, the temporal profile of this outburst is unlike giant 
outbursts in other Be/X-ray binaries that typically consist 
of either a single large intensity peak or a large peak followed 
by smaller peaks modulated at the orbital period. The initial 
large peak is usually much brighter and of much longer dura- 
tion that any succeeding peaks. XTE J011 1.2-7317 instead 
has two peaks that are similar in intensity with the first 
peak lasting about 25 days and the second lasting about 60 
days, followed by a declining tail lasting about 30 days. The 
separation between the two peaks is about 30 days. This 
double peaked structure is most likely not due to orbital 
effects because the frequency derivative is strongly corre- 
lated with the pulsed flux and also shows the double peaked 
structure. Such a correlation is predicted by simple accre- 
tion theory for systems accreting from a disk. If the primary 
source of spin-up was orbital, we would not expect a strong 
correlation with the pulsed flux and we would expect to see 
both increases and decreases in pulse frequency as the pulsar 
moved relative to the observer. 

3.2 The nebulosity 

The low surface brightness of the nebulosity, and the pres- 
ence of faint field stars superimposed upon it and nearby, 
combine to make a morphological interpretation difficult. 
The Ha image of XTE J011 1.2- 7317 taken from the ESO 
data shows the nebulosity clearly but it is unclear whether 
some of the brighter patches are stars or enhanced nebular 
emission. However, because the Ha filter bandpass is effec- 
tively a subset of the broader R filter, the stellar continuum 
emission present in the Ha images can be largely removed 
using an R band image. 
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for XTE J0111.2-7317. 
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Figure 7. Ratio of the 2-60 kcV pulsed flux in mCrab units 
measured with the RXTE PCA to the 20-50 kcV pulsed flux in 
mCrab units measured with BATSE. 

The final result shows that the emission is smoother 
than would appear from the Ha image alone, though two 
slightly brighter features are weakly visible to the SSW. 

3.2.1 Image characteristics 

Initial analyses of the nebula used CCD imaging from the 
SAAO 1.0m provided the 'first look'. Though images were 
obtained in all filters, only the Ha image shows nebulosity. 
Though H/3 emission is present (Table 2), it is apparently 
too weak to register in the image, due both to the large 
amount of continuum also transmitted and to the poor re- 
sponse of the CCD in the blue. 

However these data were to some extent superseded by 
the archival images from the ESO 2.2m telescope due to 
their superior depth and resolution. Images in Ha, Ha con- 
tinuum, [OIII] and [SII] were dearchived, but again only 
Ha clearly displayed nebulosity (Figure 8). Continua were 
removed for both the [SII] and [OIII] images. Unfortunately 
only Ha continuum images were available to model the con- 
tinuum component; this was not a problem with [SII] (at 
A6717,6731 c.f. Ha at A6563) because of the minimal spec- 
tral separation, but the [OIII] A5007 continuum was not ade- 
quately represented by the Ha continuum image, thus leav- 
ing significant stellar images and preventing the detection 
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Figure 8. Continuum subtracted Ha image of the region around 
XTE J0111. 2-7317. Image from ESO La Silla 2.2m archive, 
smoothed with a a = 1" Gaussian filter. Bright and dark lines 
are due to bad CCD columns. Image size is -~40 x 50 arcsec. 

of the small [OIII] nebula (Section 3.3.4). The FWHM of 
XTE J0111.2-7317 was measured both in the un-subtracted 
and continuum-subtracted [OIII] images and found to be 
consistent with the other field stars. [SII] appears to show 
extremely weak emission coincident with Ha. 

Thus analyses of the structure of the nebula in the light 
of these species were undertaken spectroscopically (see be- 
low). Using the maximum diameter D of 20" at a distance 
of 63.1pc one arrives at a diameter of 6.1pc for the Ha neb- 
ulosity. 

Of the possible nebular types considered the simplest 
is a conventional HII region, photoionised by XTE J0111.2- 
7317. A second possibility is that the nebula is an SNR. In 
view of the association with the BeXRB system, it would 
be extremely interesting to identify the SNR corresponding 
to the formation of the system's neutron star. This would, 
amongst other things, enable the system's age to be deter- 
mined. 

Another possibility which exists is that the nebulosity 
is a bowshock. Such structures arise when ram pressure due 
to supersonic motion of a star confines its expanding wind. 
The head or apex of the structure defines a point where the 
momentum of the expanding stellar wind balances that of 
the oncoming interstellar medium. Bowshocks are common- 
place around OB runaway stars. Indeed, according to the 
standard theory of BeXRB formation, such systems MUST 
be runaways (van den Heuvel et al, 2000). Therefore bow- 
shocks could be commonplace around BeXRBs. The nature 
of the nebula is discussed below. 

3.2.2 Nebular spectrum extraction 

No nebulosity had been observed when the SAAO spectra 
were taken so the exposures were calculated to properly ex- 
pose for the Be star. However since the slit of the spec- 
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trograph was aligned E-W subsequent examination of the 
raw spectral images revealed weak nebular emission lines off- 
set to one side of the stellar continuum. This one-sidedness 
stems from the asymmetrical morphology of the nebulosity. 

Two spectra of the source were acquired, the second be- 
cause the intended star largely missed the slit on the first. 
However both contained strong nebular features, and so were 
co-added. The exposures were sequential and as a pair were 
bracketted by arc exposures; no wavelength shift was de- 
tectable. All further reduction took place on this co-added 
image. 

The SAAO blue exposure similarly shows only very 
weak nebular H/3 emission. In view of the existence of a 
far higher S /N blue spectrum from La Silla, no further work 
was undertaken on the blue SAAO spectrum. 

Thus the STARLINK task EXTRACT was used to ex- 
tract the region immediately adjacent to the star including 
these emission features; sky subtraction was from a region 
10 pixels below. Whilst sky subtraction was used to identify 
and reject telluric and other non-nebular lines, actual mea- 
surement used subtraction of a constant background to sig- 
nificantly reduce noise. This also enabled telluric lines to be 
used as wavelength standards, for example the 3.682A sep- 
aration of telluric Ha from its SMC component provides 
a highly accurate velocity. Care was taken to exclude any 
stellar (including Be circumstellar disc) flux. Table 2 lists 
the nebular emission lines obtained from the SAAO red and 
ESO blue spectra. 

The ESO blue spectra taken on 3 November 1999 were 
extracted using the same method. 

Dividing the Ha flux density integrated along the slit 
by the approximate scale of 8" for the object one derives 
a Ha flux density of ~ 20 x W~ 16 erg/cm 2 / s/arcsec 2 , 
which is close to the value found for the bowshock as- 
sociated with Vela X-l (Kaper et al, 1997) of ~ 10 x 
10^ 16 erg /cm 2 / s/arcsec 2 . 

Because of the small number of counts from the nebular 
lines, the extracted spectra are a spatial average over all the 
nebulosity which fell upon the spectrograph slit. Without 
higher S/N data this is the only way of extracting spectra of 
usable quality. Discussion on the spatial distribution of the 
different emission lines is presented below. 

3.2.3 Line fitting 

In order to derive the most accurate line fluxes the Emis- 
sion Line Fitting (ELF) suite within DIPSO was employed. 
The FWHM values were specified from measurement of the 
H[3 line to reduce the number of free parameters and thus 
errors. The final fit is shown in Figure 9. 

The velocities are derived from the shift relative to labo- 
ratory rest wavelengths, and further corrected to a heliocen- 
tric reference frame using the STARLINK program RV. The 
positive values refer to a movement away from the observer, 
i.e. a redshift. 

3.2.4 Line flux calibration 

The presence of the stellar continuum alongside the nebu- 
lar spectrum allows an absolute flux calibration to be per- 
formed. Such calibration is important for establishing reli- 
able line ratios for diagnostic techniques. 



[Oil] line ratio 



1 

_n 1-1 




III 
^nn„ n n n nOh^^^ 














Residuals 


JV x — 















_ 2 I I I I I I I I I 

3690 3700 3710 3720 3730 3740 3750 3760 3770 



Wavelength 

Figure 9. Gaussian fit to the [Oil] A3727,3729 lines seen in the 
nebulosity around XTE J0111. 2-7317. 

The reddened theoretical Kurucz model photospheric 
spectrum (Kurucz, 1979) fitted by CHR was further red- 
dened to allow for the circumstellar reddening. This was 
normalized using the photometry at B and R for the blue 
and red spectra respectively. The extracted stellar spectrum 
was divided into this and fitted with a polynomial yielding 
the required scaling factor at each wavelength to convert 
from counts (the units of the line fitting) to erg/cm 2 /s/ A. 
Clearly this calibration technique is dependent upon the Ku- 
rucz atmosphere being an accurate model of the true stellar 
spectrum. Note however that only the shape of the spectrum 
within each spectral range is important; the normalization, 
which is performed using photometry approximately in the 
centre of each range, ensures that the flux density is correct 
at this point and thus near the mark over the entire spectral 
range. 

3.3 Spectral results 

3. 3. 1 Temperature. 

The O line strengths constrain the nebular electron temper- 
ature to below 25-40, 000K; in fact standard nebular theory 
provides a mechanism to maintain the temperature in the 
range 5000 - 10000K. Collisions of high energy electrons with 
ionised O, Ne, S and Fe excite these ions into the metastable 
states which subsequently decay by photon emission. This 
mechanism can be seen to be at work from the high inten- 
sity of O and S forbidden lines. This effectively limits the 
electron temperature to 10000K with a mean value of 6700K 
(Kitchen 1987). The consistency with which this mechanism 
has been seen to provide a 'thermostat' controling the tem- 
perature of nebulae provides encouragement that it also ap- 
plies in this case. 

3.3.2 Nebular Density. 

[Oil] and [SII] both emit close pairs of optical lines whose 
ratios are diagnostic of N e , the electron density (Saraph & 
Seaton, 1970). Using ratios of lines produced by the same 
species is appealing as uncertainties based on abundance, de- 
gree of ionization and the many other variables cancel out. 
This diagnostic tool is based upon the increasing effect of 
collisional de-excitation in preventing the emission of a pho- 
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Species 


A 


Rest A 


Flux per arcsec 


Velocity 


FWHM 




(A) 


(A) 


(10- 16 er g/ cm 2 /s) 


(km/s) 


(km/s) 


ESO blue 












S600-5900 "A 












[Oil] 


3728.12 ±0.14 


3726.16 


25.3 ±3.7 


145.3 ± 10.8 


176 ± 20 


[Oil] 


3730.91 ±0.14 


3728.91 


37.4 ±4.0 


148.4 ± 10.8 


176 ± 20 


H7 


4343.10 ±0.43 


4340.47 


10.0 ±2.7 


169.7 ±29.4 


217 ±60 


H/3 


4864.31 ±0.10 


4861.33 


17.7 ± 1.3 


172.0 ±6.0 


189 ± 19 


[OIII] 


4961.89 ±0.23 


4958.91 


2.8 ± 1.2 


168.1 ± 13.6 


171 ± 39 


[OIII] 


5009.82 ± 0.23 


5006.84 


7.7 ±1.4 


166.4 ± 13.5 


170 ± 39 


SAAO red 












6180-6930A 












Ha 


6566.28 ± 0.02 


6562.76 


150.5 ± 3.5 


168.3 ±8.3 


66 ±2 


[SII] 


6720.03 ± 0.06 


6716.47 


25.1 ±3.5 


163.7 ±5.7 


47 ±6 


[SII] 


6734.41 ± 0.06 


6730.85 


18.0 ±3.2 


163.4 ±5.7 


47 ±6 



Table 2. List of observed discrete emission lines from nebulosity around XTE J0111.2-7317. Wavelengths from Allen (1976), Meinel 
(1965) and Takami et al. (2001). 



ton by an excited ion in a metastable state, as one considers 
progressively higher densities. 

In the neighbourhood of the critical densities from 10 2 — 
10 4 cm -3 accurate densities may be derived, but outside of 
this range the results are indistuinguishable from either or 
oo on each side of the critical density. The results in Table 3 
(derived from line fits shown in Figure 9) place an upper 
limit of iV e < 200cm" 3 . 

3.3.3 Line velocities. 

The heliocentric velocities of the nebular lines are all consis- 
tent with a systemic SMC velocity of 166km/s, excepting the 
[Oil] A3726,3729A. However, these lines lie at the extreme 
end of the spectral range where the wavelength calibration 
is less accurate, so this is not believed to represent a real 
velocity shift. 

The FWHM are all consistent with instrumental broad- 
ening, with values comparable to those of the calibration arc 
lines. The resolution of the blue spectrum is such that struc- 
ture with a velocity dispersion of up to ~ lOOkm/s could lie 
unresolved, whereas the red spectrum limits such structure 
to ~ 30km/s. This result argues strongly against a SNR hy- 
pothesis for the nebulosity; even 'old' SNR display expansion 
velocity dispersions of several hundred km/s (Lozinskaia, 
1980). 

3.3.4 Spatial distribution of nebular lines. 

The full known extent of the nebulosity is defined by the 
narrowband Ha image seen in Figure 8. Whether this rep- 
resents the true distribution of material in the structure or 
is more a reflection of the state of ionization can best be 
inferred by looking at the distribution of emission from dif- 
ferent atoms and ionization states. 

Representative cross-sections along the spatial axis of 
the spectral images are shown in Figure 10, taken through 
the mid-point of the detected emission lines with an appro- 
priate width of pixels sampled (mostly 3) depending upon 
the FWHM. These plots thus show the spatial distribution 
of emission from each emission line; with the caveat that the 
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Figure 10. Cross sections through the full widths of the [Oil] 
and [OIII] emission lines. The stellar contribution has been sub- 
tracted, and is shown for reference. Profiles have been Gaussian 
smoothed with a = 0.7. 



Be star produces a peak at its spatial location. This arises 
from both photospheric emission and from the circumstellar 
disc. These peaks were removed by subtracting an extracted 
profile derived by interpolation of the continua adjacent to 
each spectral line, in a manner analogous to sky subtrac- 
tion in conventional spectral reduction. This was relatively 
successful for those lines which have no photospheric or cir- 
cumstellar counterpart: [SII], [O II] and [OIII], though [O II] 
was affected by the rapid non-linearly declining response of 
the CCD at short wavelengths. 

In the case of the Balmer lines the combined spatial 
profile from the photosphere and circumstellar disc has the 
same FWHM in the spatial direction as the continuum, but 
must be scaled to account for the circumstellar emission. 

The orientation of the slit for the ESO blue spectrum 
- [OIII], [Oil], H/3 and H-y data is NW - SE through the 
star and the brightest part of the nebulosity, thus in the 
context of the bowshock scenario it is usefully aligned along 
the direction of the standoff distance I. The [OIII] line was 
expected as it is observed to be amongst the strongest lines 
in many nebulae, including other bowshock nebulae, HII re- 
gions and PNe. 
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Narrow band imaging of the field has not produced a 
detection at [OIII] A4959,5007 and in the spectral data it is 
detected only weakly from the nebular regions brightest in 
H[3. It is however observed as a strong feature from regions 
much closer to the star itself, to the extent that extraction 
of the stellar spectrum shows strong features at the afore- 
mentioned wavelengths (Figure 2), whereas extraction of 
the neighbouring nebular spectrum detects them only very 
weakly. This can be understood with reference to Figure 10 
which shows spatial profiles of the [Oil] and [OIII] emission 
lines along the spatial axis. 

Whilst the [OIII] emission peaks close to the stellar lo- 
cation and rapidly drops off, H/3 is observed to considerably 
larger distances. The [OIII] emission profile is much more 
symmetric around the Be star, with only minor enhancement 
in the direction of strongest Balmer emission. Additionally 
[OIII] is seen from the NW side of the Be star where there is 
negligible Balmer emission. Where Hp and H'y are brightest 
there is scarcely any [OIII] emission at all. 

The HP, H'y and [Oil] emission appear to be coinci- 
dent, though few conclusions can be drawn from the weak 
Hy detection. All are seen only to the SE side of the Be 
star, in a manner totally consistent with the distribution 
seen from the Ha image in Figure 8. Detection from pixels 
75.5 to 92.5 implies an angular diameter of 13.9", or 4.3 pc. 

The SAAO red spectrum - Ha and [SII] - shows the 
distribution of nebular Ha emission. The size of the nebu- 
losity appears to be the same as that of Hp and Hy, with 
detection out to slightly larger radii because of the much 
stronger signal. [SII] appears to trace the Ha exactly. 

3.4 Nebular classification. 

If a spectral classification of B0.5-B1 (CHR and Covino et 
al, 2001) is accepted it becomes possible to calculate the 
Stromgren radius within which a uniform cloud of H is com- 
pletely photoionized (note this does not apply to the hol- 
low bowshock model) and would be expected to emit Ha. 
This radius is sensitive to spectral class around BO, in that 
this classification sensitively determines the Lyman contin- 
uum flux Qo depending upon the details of the stellar model 
used. Vacca et al. (1996) tabulate Qo for O3-B0.5V stars, so 
though one must perform an extrapolation to yield a re- 
sult for Bl, the function is a smooth one and so the result 
expected to be accurate. Incorporating the Martins et al. 
(2002) downward revision by 40% of Vacca et al.'s Qo, one 
arrives at 4.8 x 10 47 at B0.5V and 2.4 x 10 47 at B1V. 

Suitably modifying the derivation of equation 6A.1.10 
from Kitchin (1987), we arrive at the formula 

Rs = [ ; — : fj-oF ) 3 (metres) (2) 

\4ttN? x 6.5 x 10- 15 T e -°- 8B y 

Using the observed nebular radius of 3.05pc requires 
iV e = 7.6cm" 3 at B0.5 or only iV e = 5.1cm" 3 at B1V, at the 
lower end of observed nebular densities, and in agreement 
with the results of Section 3.3.2 using forbidden line ratios. 
The use of SMC metallicities, reducing line blanketting ef- 
fects, raises Qo and thus requires slightly higher densities, 
but this effect is only of the order of a few %. 

This suggests the scenario that the nebulosity is simply 
an HII region, perhaps part of a cloud adjacent to and ion- 



[OII] 3729:3726 1.48 ± 0.27 

[SII] 6716:6731 1.40 ± 0.31 

[SII] :Ha 0.286 ±0.031 

H/3: [OIII] 1.69 ±0.32 



Table 3. Diagnostic line ratios from nebulosity around XTE 
J0111. 2-7317. 

ized by XTE J0111.2-7317. If so, and considering the fact 
that the calculated densities are very much at the lower end 
of the range for HII regions, why is this structure detectable 
at all? Most OB stars rid their neighbourhood of gas early 
in their existence via powerful winds. The answer may lie 
in the fact that as a BeXRB, XTE J0111. 2-7317 has a large 
space motion (van den Heuvel et al, 2000) and has im- 
pinged upon a cloud which it has subsequently ionized. If 
this cloud increases in density towards the south-east (lower 
left in the images), the slightly smaller radius in this direc- 
tion can be explained with the Stromgren sphere argument 
(Rs = fn(Nij)), or as an ionization front propogating into 
the cloud. 

3.4-.1 The bowshock theory 

Stellar bow-shocks form from the ram-pressure interaction of 
the local interstellar medium with the stellar wind; when the 
star's velocity exceeds the sound velocity in the ISM material 
is swept up into a shock front which is detected primarily 
using one of two techniques: narrowband imaging (mostly 
Ha) and IR imaging, mostly using IRAS data. Whereas 
Ha imaging takes advantage of excited H, the IRAS flux 
(usually most pronounced in the 60^im band) arises from 
thermal radiation from swept-up dust originating from the 
star, heated by shock interactions. 

The appearance of a bowshock is clearly heavily depen- 
dant upon the angle that the star's velocity vector makes 
with respect to the plane of the sky. Classic parabolic struc- 
tures such as Vela X-l and a Cam are seen when the star's 
motion lies almost entirely within this plane. Less clearly de- 
fined structures are more commonly seen, though enhanced 
brightness on one side, such as is observed in XTE J0111.2- 
7317, is usually apparent (Noriega-Crespo et al, 1997). 

Thus the morphology of a bowshock is such that the 
head or apex of the structure is well defined if viewed side-on, 
enabling a 'standoff distance' I to be determined; there is no 
sharply defined trailing edge. The size of the structure is thus 
easily characterised only for side-on specimens. In the case 
of XTE J0111.2-7317, the standoff distance is approximately 
3pc, within the range found for other bowshocks tabulated 
in Table 4. 

Comparing the 60/im IRAS image of the a Cam bow- 
shock in with the XTE J0111.2-7317 nebula reveals an ex- 
tremely similar appearance. Both nebulae lie within an ap- 
proximately parabolic perimeter, and though a Cam itself is 
not visible in the IRAS image, both stars lie slightly inside 
the radius of curvature of the nebulosity, agreeing with the 
theoretical radius of curvature for a stellar wind bowshock 
of (5/3)/. 

Without doubt a measured space velocity for XTE 
J0111.2-7317 would be strong evidence one way or the other. 
The radial velocity has been investigated by cross correlat- 
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Object 


Spectral class 


Standoff dist.(pc) 


Space vel.(km/s) 


XTE J0111.2-7317 


BOVe 1 


3.0 




Vela X-l 


BO.blbe 2 


0.48 


~ 50fcm/s 8 


a Cam 


09.5Iae 3 


5.1 6 


48 7 


0623+71 


CataclysmicVar. 4 


0.08 


~ lOOkm/s 


Betelgeuse 


M2Iab 5 


0.8 


56 



Table 4. Sizes of some stellar bowshocks and spectral classification of their associated stars. Sources l:This work, 2:Kaper et al. (1997), 
3:Noriega-Crespo et al. (1997), 4:Hollis et al. (1992), 5:Noriega-Crespo ct al. (1997), 6:van Buren and McCray (1988), 7:Stonc (1979), 
8:Comeron & Kaper (1998) 



ing the blue ESO spectrum with that of a velocity standard, 
but errors of 50km/s and a result consistent with zero were 
sufficiently large to preclude any conclusions being drawn, 
probably because of the low S/N of the spectrum. Studies 
of specific lines would seem to suggest a limit on v r of less 
than lOkm/s. The bowshock appears however to be nearly 
edge on, in which case most of the stellar velocity is oriented 
in the plane of the sky and a modest or zero v r would be 
expected. 

Comeron & Kaper (1998) have performed numerical 
simulations of bowshocks in a number of cases. One of their 
findings is that in low velocity cases (v* <C lOOkm/s), of 
which XTE J0111. 2-7317 is an example as a BeXRB, the 
bowshock becomes much thicker as the assumption of in- 
stantaneous cooling, implicit in the above equation, fails. 
The resulting thick layer of hot gas, bounded on one side 
by the stellar wind, and on the other by the ISM, certainly 
can explain the broader appearance of XTE J0111.2-7317's 
nebula compared to that of Vela X-l. 

A second consequence is that peak density occurs at 
1.5 — 21, somewhat away than the balance point I predicted 
by simple momentum balance. Rearranging equation 16 of 
Raga et al. (1997) to incorporate metallicity dependances 
(Vink et al, 2001) of M and vw gives the following equation 
balancing the ram pressure of the ISM with that of the ex- 
panding stellar wind at the standoff point, a distance I from 
the star. 



N 



ISM 



M Z=1 Z 0S2 



v w 



A-wPmnv 2 



(3) 



Using a solar metallicity mass loss rate of 3 x 
lO^MOj/r" 1 (Vink et al, 2000), a wind velocity of 
lOOOkm/s, the assumption of a pure hydrogen ISM, system 
velocity V w = 15km/ s (mean BeXRB velocity), and I — 
3.0pc/(1.5 - 2) requires an ISM density of 0.02 - 0.034cm" 3 . 
This value is low for any ISM, particularly in a gas rich 
galaxy like the SMC where densities of less than 0.1cm" 3 
are uncommon. The assumption of an unusually high wind 
speed and/or mass loss rate can reconcile this situation, but 
both seem somewhat unlikely. 

Though this is clearly some way from an exact method 
of calculating the ISM density, it seems safe to conclude 
that it must be low if we believe the bowshock interpreta- 
tion. We are therefore led to questioning why this is the 
only one known around a BeXRB. The requirement of su- 
personic velocity may be the reason; the typical ISM sound 
speed of vsctN 1 ^ 2 ~ 10/cm/s is only slightly below the mean 
BeXRB velocity of 15km/s (van den Heuvel et al, 2000), so 
perversely it could be only the BeXRB systems travelling 



through relatively tenuous regions that are capable of pro- 
ducing bowshocks. 

The argument in Section 3.4 concerning the Stromgren 
radius around the B star does not apply to such thin shells. 
Equation 21 of Comeron and Kaper (1998) considers the flow 
of material into the bowshock to find the required ionizing 
flux - using the values given above, the flux necessary to fully 
ionize the material in the shock front is only 4 x 10 42 sec _1 , 
which is several orders of magnitude less than a B0.5-B1 
star provides. Thus we expect to see Ha emission from the 
bowshock. 

Whilst the [Sll]:Ha ratio of 0.286 is below the threshold 
of 0.4 usually accepted as clear evidence of shock excitation, 
it is significantly above the typical values for photoionized 
nebulae which mostly lie below 0.2. This suggests the possi- 
bility that a combination of photoionization and shock ex- 
citation may be at work. 



3.4-2 Comparison with other bowshock spectra 

The discovery and study of bow-shocks around high velocity 
stars is a relatively new facet of astronomy. Few published 
spectra exist - spectra of only two other stellar wind bow- 
shocks have been found in the literature: Vela X-l (Kaper 
et al, 1997) and a CV system 0623+71 (Hollis et al, 1992). 
Both sources show strong Balmer emission as well as [OIII] 
and [SII] lines, in common with XTE J0111.2-7317. Addi- 
tionally the 0623+71 spectra show strong [Oil] lines much 
like XTE J011 1.2- 7317. The only significant difference is 
the absence of [Nil] lines in the XTE J011 1.2- 7317 spec- 
tra, though the Vela X-l spectra show them so weakly that 
if present in XTE JO 11 1.2- 73 17 at a similar level they would 
not have been detected. Thus the Vela X-l nebular spectrum 
is not significantly different to that of XTE J011 1.2- 7317, 
strengthening the bowshock hypothesis. 

Essentially the XTE J011 1.2- 7317 nebular spectra are 
qualitatively similar to both of the bowshock spectra shown, 
in particular that of Vela X-l. This is to be expected as the 
XTE J0111. 2-7317 system far more closely resembles Vela 
X-l than 0623+71. 

The presence of [OIII] far from the presumed shock in- 
teraction suggests photoionization as the mechanism for its 
excitation rather than shock proceses; Noriega-Crespo and 
Garnavich (1992), in their studies of bowshocks of Herbig- 
Haro objects, confirm earlier findings (Dopita et al, 1984) 
that shock velocities of 80km/s are required to provide suf- 
ficient excitation to produce [OIII] A4959,5007, and spectra 
do not reveal such velocities to be present. Bowshock theory 
also does not predict this order of velocity. 
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3.4-3 The HII region hypothesis. 

A high [Olll]:H/3 ratio is a diagnostic of photoionization 
from the Be star's UV flux, comparable to the strong [OIII] 
lines seen in planetary nebulae caused by the white dwarf's 
UV flux. The observed ratio (Table 3) is somewhat higher 
than seen in the two available bowshock spectra and is typ- 
ical of HII nebulae and at the large end of the distribution 
for SNR (Fesen et al, 1985; Blair & Long, 1997). 

The observation of a small [OIII] nebula coincident with 
the central star raises a possible problem for the bowshock 
theory. If this material truly lies within the larger Ha neb- 
ula, it would appear to argue against the bowshock theory, 
as the interior of its paraboloidal shell is considered to be 
empty, excepting the expanding rarified stellar wind. The 
simulation work of Comeron & Kaper (1998) allows for some 
material to pervade this region, but probably not as close to 
the star as we observe. Therefore in this scenario the [OIII] 
emission must lie on the surface of this shell. In the case 
of XTE J0111. 2-7317 however, the chances of one patch of 
emission lying so well centred on our line of sight to the 
star seem remote - the best explanation for this must be 
that the [OIII] emission arises physically within the bounds 
of the Ha nebula, particularly as a Stromgren sphere argu- 
ment can account for the relative sizes and locations of these 
features in this scenario. 

Section 3.3.4 shows that the HII, [SB] and [Oil] nebu- 
lae are approximately coincident. Whereas this flows natu- 
rally from a bowshock theory (one is simply measuring the 
standoff distance I), it can also be understood in terms of 
Stromgren spheres. 

In view of the abundance of nebulae within 30' of XTE 
J0111.2-7317, the association with a gas cloud does not cause 
any problems. SHASSA (Southern B-Alpha Sky Survey At- 
las) (Gaustad et al, 2001) shows large BII regions close by 
and diffuse emission to be extremely pervasive. 

The high space velocity which all BeXRB possess could 
simply have enabled the star to drift into an interstellar 
cloud and ionize it; it could not have been in this position 
long as the surrounding gas would have been blown away 
by the stellar wind. The mean space velocity of BeXRBs of 
15km/s (van den Heuvel et al, 2000) corresponds to lpc in 
11,000 years. 

3.44 SNR hypothesis 

Obviously, the presence of the NS in a BeXRB requires there 
to have been a SN explosion at some stage in the system's 
evolution. Thus the discovery of a nebula surrounding the 
system naturally raises the possibility that this may be a 
SNR, thereby enabling the age of the system in its BEXRB 
state to be reliably estimated and providing verification of 
the BeXRB formation mechanism. However, spectral evi- 
dence shows that this is not the case. 

Even a qualitative comparison of the nebular spectrum 
reveals that it is much more closely allied both to published 
bowshock (discussed above) and HII region spectra than to 
a those of SNR. In addition, several quantitative properties 
argue powerfully against it being a SNR. 

The [SII]:-ffa ratio is a good distinguishing diagnostic 
between HII regions and SNRs, and is the standard diag- 
nostic used for classifying nebulae in neighbouring galaxies 



where morphology cannot be used (Blair & Long, 1997). 
With the exception of a handful of extremely low surface 
brightness HII regions similar to and including the diffuse 
interstellar gas, the vast majority of such objects have a 
[SB] :Ha ratio of less than 0.4, and usually below 0.2. SNR 
have values normally significantly in excess of 0.5. Thus with 
a ratio of 0.286 XTE J0111. 2-7317 lies at the upper end of 
the distribution of HII regions, suggesting the possibility of 
a shock heating component in addition to photoionization 
as the dominant mechanism. 

The temperature sensitive line ratio [OIII] A5007/[OII] 
A3726 + A3729 appears to be 0.12 from Table 2, implying 
a temperature in excess of 70,000K (Figure 7 of Dopita, 
1997) this is strongly affected by the photoionized compo- 
nent at small distances from the star; no [OIII] emission is 
detected from the shock-front itself implying temperatures 
below those encountered in SNR. Similarly, whilst the [OIII] 
A4363/A5007 ratio is traditionally used to provide temper- 
atures in SNR, upper limits to the A4363 line (3 x 10~ 17 
erg/s) combined with the [OIII] line strength can only sug- 
gest temperatures below 25-40, 000K, in the realm more typ- 
ical of planetary nebulae than SNR. The negligible presence 
of ++ at the shock location in itself betrays temperatures 
of below 40,000K (Dopita 1977). Additionally the lack of 
[O I] A6300,6364 places the system well away from the realm 
of SNR in the diagnostic diagrams of Fesen et al. (1985). 

3.5 Other HMXBs 

A total of 15 other HMXBs in the Magellanic Clouds were 
examined for evidence of nebulosity, but none was found. 
Anything down to a scale of ~ 1 pc and Ha flux density 
- 20% of the XTE J0111.2-7317 nebula (which has almost 
exactly the same Ha surface brightness as the Vela X-l bow- 
shock) would have been detected. In view of the consistency 
of the surface brightness of these structures and typical sizes 
it seems probable that any other similar structures would 
have been found. Thus the frequency of such objects in our 
sample is 1 in 15, or 6.7%. This statistic is in exact agreement 
with the results of Huthoff and Kaper (2001) who searched 
for IR bowshocks around 15 galactic HMXRBs, and found 
none except for the well-known case of Vela X-l (Kaper et 
al, 1997). 



4 CONCLUSIONS 

The observed characteristics of the nebulosity surrounding 
XTE J0111. 2-7317 have been examined in the context of 
three possible nebular types: SNR, bowshock and HII re- 
gion. Of these, the SNR hypothesis is excluded by several 
line ratios and velocity dispersion. Such properties provide 
no discrimination between bowshocks and HII regions, but 
the spatial distribution of line emission from different species 
matches predictions for an HII region to a high degree, whilst 
providing only a poor match to the expectations for a bow- 
shock. Overall these results provide good evidence that the 
nebulosity surrounding XTE J0111.2-7317 is a conventional 
HII region. 

The x-ray outburst from the 31-s X-ray pulsar in XTE 
J0111.2 observed with BATSE and RXTE in 1998-1999 ex- 
hibited characteristics typical of a giant or type II outburst 
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in a BeXRB including large spin-up rates, Lx > 10 erg/sq- 
cm-s, and a correlation between spin-up rate and pulsed flux. 
However, the temporal profile of the outburst was unusual, 
consisting of two similar intensity peaks, with the first peak 
of shorter duration than the second. 

Thus, while apparently unusual, XTE J011 1.2- 7317 ap- 
pears, in fact, to be a normal SMC BeXRB embedded in a 
locally enhanced ISM which it has photoionised to create an 
HII region. 
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